The QR iteration method for Hermitian quasiseparable matrices of an arbitrary order
نویسندگان
چکیده
The QR iteration method for tridiagonal matrices is in the heart of one classical method to solve the general eigenvalue problem. In this paper we consider the more general class of quasiseparable matrices that includes not only tridiagonal but also companion, comrade, unitary Hessenberg and semiseparble matrices. A fast QR iteration method exploiting the Hermitian quasiseparable structure (and thus generalizing the classical tridiagonal scheme) is presented. The algorithm is based on an earlier work [6], and it applies to the general case of Hermitian quasiseparable matrices of an arbitrary order. The algorithm operates on generators (i.e., a linear set of parameters defining the quasiseparable matrix), and the storage and the cost of one iteration are only linear. The results of some numerical experiments are presented. An application of this method to solve the general eigenvalue problem via quasiseparable matrices will be analyzed separately elsewhere. AMS classification: 15A18; 15A57; 05E35; 65F15
منابع مشابه
Eigenstructure of Order-One-Quasiseparable Matrices. Three-term and Two-term Recurrence Relations
This paper presents explicit formulas and algorithms to compute the eigenvalues and eigenvectors of order-one-quasiseparable matrices. Various recursive relations for characteristic polynomials of their principal submatrices are derived. The cost of evaluating the characteristic polynomial of an N × N matrix and its derivative is only O(N). This leads immediately to several versions of a fast q...
متن کاملAn iterative method for the Hermitian-generalized Hamiltonian solutions to the inverse problem AX=B with a submatrix constraint
In this paper, an iterative method is proposed for solving the matrix inverse problem $AX=B$ for Hermitian-generalized Hamiltonian matrices with a submatrix constraint. By this iterative method, for any initial matrix $A_0$, a solution $A^*$ can be obtained in finite iteration steps in the absence of roundoff errors, and the solution with least norm can be obtained by choosing a special kind of...
متن کاملStability of QR-based fast system solvers for a subclass of quasiseparable rank one matrices
The development of fast algorithms to perform computations with quasiseparable matrices has received a lot of attention in the last decade. Many different algorithms have been presented by several research groups all over the world. Despite this intense activity, to the best of our knowledge, there is no rounding error analysis published for these fast algorithms. In this paper, we present erro...
متن کاملComputing the Matrix Geometric Mean of Two HPD Matrices: A Stable Iterative Method
A new iteration scheme for computing the sign of a matrix which has no pure imaginary eigenvalues is presented. Then, by applying a well-known identity in matrix functions theory, an algorithm for computing the geometric mean of two Hermitian positive definite matrices is constructed. Moreover, another efficient algorithm for this purpose is derived free from the computation of principal matrix...
متن کاملA Study of the Hr and Extended Hr Methods for the Standard Eigenvalue Problem
The QR method is a very eecient method for computing the spectrum of Hermitian tridiagonal matrices since the tridiagonal form is preserved over the iterations. For non-Hermitian tridiagonal matrices the QR method destroys the tridiagonal form. In this report we study two methods, the HR and the XHR methods, that preserve tridiagonal form for pseudo-Hermitian matrices. We also report results fr...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2005